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Some methods for the numerical computation of two-loop non-infrared 
vertices are reviewed. A new method is also proposed and compared to the 
old ones. Finally, some preliminary results are presented, concerning the 
evaluation of the fermionic corrections to sin 2 throught the described 
techniques. 



1. Introduction 

In the forthcoming experiments, the validity of the Standard Model 
will be tested with high precision. In addition to the direct search of the 
Higgs Boson, important quantities will be measured in the future colliders, 
providing a good test of the Model. This of course pushes the theorists to 
compute these observables at the same degree of precision. For example, 
the mass of the W boson, whose present value is Mw = 80.426 ±0.034 GeV 
(U), will be measured with an expected error of 15 MeV at LHC and 6 
MeV at the ILC. Or the effective leptonic weak mixing angle (sin 2 9^ = 
0. 23150(16), P]) will be known with an absolute precision of 10 -5 at the ILC. 
To get a similar theoratical incertenty, we have to improve the calculation 
in perturbation theory beyond the one-loop level. The computation of two- 
loop Feynman diagrams is a hard task. The pure analytical techniques are 
very efficient when few mass scales are present (see for example j5] or [3]), 
but seem to be unable to deal with the complete set of two-loop diagrams 
in the Standard Model (where more scales come into the game). For this 
reason we were led to abandon the analytical way in favor of a numerical 
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evaluation of multi-loop diagrams. The goal of the numerical approach is 
to express any diagram in terms of smooth integrals. 



2. Standard BT relation 



The Bernstein-Tkachov theorem jl] tells us that for any finite set of poly- 
nomials Vi(x), where x = (xi, . . . vector of Feynman parameters, 
there exists an identity of the following type: 



T(x,d)l[Vt + \x)=Bl[vnx), 



(1) 



where T is a polynomial of x and di = d/dxi\ B and all coefficients of T 
are polynomials of fi{ and of the coefficients of Vi(x). If the polynomial V is 
of second degree we have a master formula, again due to F. V. Tkachov [I]. 
We write the quadratic V as: 



V(x) = x t H x + 2K t x + L, 



(2) 



where x = (x±, ...,x n ), H is an n x n matrix, K is an n vector. The solution 
to the problem of determining the polynomial T is as follows: 



M + l' 



V 



X 



B = L-K l R- l K 



X 



-H- 1 K. 



Therefore we have: 



V(x) = i 



1 + 



This is the standard BT relation for quadratics. 



V» +l {x). 



(3) 
(4) 

(5) 



3. Strategy 

The standard BT relation is very usefull for the computation of one- 
loop diagrams and also some two-loop configurations (see [3], 0) 0) |3|> 
[Hj). The strategy is the following. First of all, the diagram with N external 
legs and a tensorial structure of rank n is decomposed in form factors: 

= J2 G N (T) T w '-'^. (6) 

T 

Here T^ 1 '"''^™ are all possible tensors of rank n that can be obtained by 
combining the N — 1 independent external momenta of the diagram and the 
metric 8 ull . 
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For one- loop diagrams, if we write each form factor in the parametric 
space, we always obtain a result of this form (see [7j for details): 

/ 2 \e/2 

Gjv= m F ( n + fj/^-iO*) P(x)V{x)- n ^/ 2 nGiV, (7) 



where 

dS m (z) = dzi / dx 2 ... / cfe m , (8) 
jo jo jo 

/x is the mass scale and P and V are polinomials in the Feynman param- 
eters x = (x\, . . . ,Xjv_:i). In particular V(a;) is always a quadratic of the 
type Eq.pj). The T function contains the UV pole (if present) of the dia- 
gram. Its argument is always equal to the exponent of V with opposite sign 
and therefore, for every UV divergent one-loop diagram, V has a vanishing 
power. 

The goal is to express G N in terms of smooth integrals to be integrated 
numerically. For n = (which corresponds to UV divergent form factors), 
we can simply perform a Laurent expansion around e = to get just smooth 
integrands of the type: 

P{x)\n k V{x). (9) 

For n > 1 the idea is to apply the BT relations Eq.© to "raise" the power 
of V of one unit and then integrate by parts to get rid of the derivatives. 
Then the procedure is repeated for all integrals that are generated, until 
the power of V becomes — e/2, and, at this point, we proceed as in the case 
n = 0. 

The procedure is clear and leads to smooth integrals at the price of 
introducing the denominator B which of course can vanish somewhere in 
the phase space. It can be prooved that the zeros of B correspond to the 
leading Landau singularity of the diagram, but the singular behaviour is 
usually overestimated by the BT procedure (see Q)- This means that in 
the region B ~ all terms generated by Eq.© cancel one another, giving 
rise to numerical instabilities. For this reason it would be good to find a 
new relation, which should contain the real divergent behaviour for B = 
and therefore should remain stable also for B ~ 0. 

In the two-loop case the form factors are classified counting the number 
of propagator of each loop: 

G N (T) G abc (T), (10) 

where a is the number of propagators (with momentum qi) which belong 
just to the first loop (the one with the smallest number of propagators), b 
is the number of propagators (with momentum 52) which belong just to the 
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second loop and c is the number of propagators (with momentum q\ — q 2 ) 
which belong to both loops. 

Then we first parametrise the loop with momentum qi, obtaining a new 
propagator in q2 with a non integer power. The mass and the momentum 
of this new propagator depend in general on the Feynam parameters of the 
first loop. After the parametrisation of the second loop, each form factor 
takes the form ( see HQ. for details): 

G abc = A r(n + e) JdSa+^ix) J dS b {x) 

[x 2 (1 - x 2 )]- m -^ 2 yZ- 1+t,2 P(x, y) V x (y)- n - € , (11) 

where m, n G IN, fi is again the mass scale and A e is a constant regular in 
e = 0. P is a polinomial in all Feynman parameters x = (x\, . . . ,x a+c -i) 
and y = (yi, . . . ,y b ). V x (y) is a quadratic of the type Eq.Q in y, where 
now the coefficients H, K, L depend on x and have the following form: 

CW (12) 



x\ (1 — X2) k 



where h,k € IN and C is quadratic in x. The T function contains the 
overall UV pole (if present), while sub-divergencies are contained in [x 2 (1 — 
X2)]~ m ~ e l 2 y™ 1+lf//2 . Since at least one of the coefficients of V x have in 
the denominator the product x 2 (1 — x 2 ), the divergent behaviuor of the 
integrand at x 2 = and x 2 = 1 is present (generating the UV pole) just for 
m > n. On the other hand the UV divergency is generated by the behaviour 
in y b = only for m = 0. 

So, apart the UV pole coming from sub- loops, any two-loop diagram 
is the integral of a one-loop diagram whose masses and momenta depend 
on the integration variables. If we would be able to express any one- loop 
diagram in terms of smooth integrals not only with respect to the integration 
variables, but also with respect to their external masses and momenta, the 
numerical evaluation of two-loop diagrams would be a trivial task. The BT 
relation Eq.© is in general not good for this purpose, because it introduces 
the denominator B. In fact, since the coefficients H, K and LofV x fEa. (|ll|0 
depend on the Feynman variables x, the same happens for B which therefore 
generates singularities inside the x integration contour. As a consequence, 
apart some special cases where it is possible to have a factor B independent 
of any Feynman variable, the standard BT method can not be applied for 
two-loop diagrams. This is another reason to search for a new BT-like 
relation. 
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4. The new BT-like relation 

To obtain new relations it is usefull to write the quadratic V(x) in the 
following way: 

V[x) = x t Hx + 2K t x + L 

= (x t -X t )H(x-X) + B = Q(x) + B. (13) 

This formula, which defines the quadratic Q(x), can be trivially verified 
using the definition of B and X in Eq.(|4j). The basic relation satisfied by 
Q(x) is the following: 

V t d x Q(x) = -Q(x), V = -^lA. (14) 

At this point we introduce a new variable y and a new polynomial W(x, y) 
defined as follows: 

W(x,y) = Q(x)y + B, (15) 
and satisfing the following relation: 

(V t d x + yd y )W^x,y) = 0. (16) 

Next we consider the following integral 

Ip= fdyyP^Wfay), (3 > (17) 

Jo 

and compute: 

V* X Ip = [dy yP- 1 V 1 d x W»(x, y) = - [ dy / d y W»(x, y) 
Jo Jo 

= -W> M (x,l)+p [dyyP- 1 W fi (x,y) = -V fl (x) + pis. (18) 
Jo 

Using the definition of the hypergeometric function (see ^J) to evaluate 
1/3, we finally get: 

V»{x) = B» ( 1 - ^ ^ 2 F 1 (- t i , ; p + 1 ; -|). (19) 

This formula has a general validity and does not depend to the actual ex- 
pression for Q, B and V. The only relations which they must satisfy are: 

V(x) = Q(x) + B, V t d x Q(x) = -Q(x), > 0. (20) 



G 



The usefullness of this relation is evident if we consider the case \i = —1—ae. 
In this case the better choice for the free parameter /3 is f3 = 1. Using the 
expansion of the hyper geometric function around e = 

1 00 (- n A n 

2 F 1 (l + ae,l;2;x) = -- £ f— ln« +1 (l - x), (21) 

X n =0 ^ + '■ 

we obtain 

- g (1 - ^ *) ^ h,- (l + §) . (22) 

In this relation we have obtained our goal to avoid the appearence of the 
factor B in the denominator. Here, the only denominator is Q(x) which can 
vanish inside the integration contour; however its zeros are compensated by 
the logarithm, whose argument goes to 1 when Q(x) goes to 0. 

An example of the usefullness of the new relation is the evaluation of 
one-loop three-point functions. In the scalar case we have: 

°o= (y) f ^^ + ^)J/xiJ o Xl dx 2 V(x 1 ,x 2 y 1 -^ 2 . (23) 
If we insert Eq.(|22|l and integrate by parts, we simply obtain: 

Co = Ef A» tA-^"( 1 + ^ir 1 ) + °< £ >- < 24 > 

^ 2 Jo Q[t](x) V B J 

where 

Q[0](x) =Q(l,x), Q[l](x) = Q(x,x), Q[2](x) = Q(x, 0), (25) 

a = l-X 1 , a 1 =X 1 -X 2 , a 2 = X 2 . (26) 

This result for Co (which can be easily generalised for tensor integrals) 
can be also used to compute two-loop diagrams which can be expressed as 
an integral of a one-loop three-point function. We see for example what 
happens in two families of two-loop vertices. 



5. The two- loop vertex V 131 

All two-loop vertices can be classified according to six families. Their 
list is given in the appendix. Taking into consideration the U 131 vertex, 
after Feynman parametrisation it takes the form: 



V 131 = ~ (^j r (1 + e)£dx JdS 3 (y, z u z 2 ) [ x (l- x )]-^(l- y y/^ 



(27) 
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u 



-X f ) H ( z -X) + (ml-m 2 3 ) (l-y)+B, 



mt 



2 

— +, 
x 1 



rrin 



-, (28) 



where the masses are defined in the figure for V 131 in the appendix. The 
coefficient H, X and B are those appearing in the polynomial of a one-loop 
three-point funcion, with external momenta —P, p±, p 2 and masses 7713, 7714, 
7715. If we introduce Q(y, z\, z%) for the polynomial U defined by: 



U{y,z 1 ,z 2 ) = Q{y,z 1 ,z 2 ) + B 
we see that Q satisfies the following basic relation: 

y - x f ) a, 



Q = -Q. 



(29) 



(30) 



From this formula we obtain the standard BT relation and the new one 
(choosing (3=1): 



U~ 



-l-e 



1 

B 



1 



(l - y) d y y - X 1 ) d z 



u 



-l-e 



E 

n=l 



(-^) 



n-1 



l-(l-y)d y + 



2e 

y - 9 2 



(31) 



ln n 1 + 



73 



(32) 

From these equations, we see that V 131 is exactly one of those particular 
two-loop configurations for which can be found a factor B independent 
from any Feynman parameter. If this is crucial for the application of the 
standard BT method (Ea.([31|0. this would not be strictly required for the 
new method fEa.(|32[)). In addition to that the new relation appear to have 
a better behaviour near the zeros of B. 

However, for this type of diagrams, where the polynomial U is linear in 
one of the Feynman variables (y), another procedure is available. After the 
transformation y — > 1 — (1 — z\) y, we have to compute: 



dyy^iay + b)- 1 '* 
1 



2 b' 



-l-e 



2^(1 + 6,1; 1 + 1; -£), (33) 



a = (ml-ml)(l-z 1 ), b = (z t - X 1 ) H (z - X) + B. (34) 

By applying the properties of the hypergeometric function (see [5] for de- 
tails) and expanding in e we simply obtain: 



V 



131 



dz\ 



21 1 
dz 2 — In [ 1 + 



V 



(ml 



m 



+ h / ln^. C Q (P ,p 1 ,p 2 ,m 3 ,m 4: ,m 5 ), 



(35) 
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where 

V = {z t - X 1 ) H(z-X) + B, A x = -ml x{l-x) + m\{l-x) + m\ x. 

Note that V is exactly the polynomial of the Co function appearing in 
Eq. (|35|) . The numerical results for the three methods are then compared in 

Fig.m 
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Fig. 1. Numerical results for the V 131 family: V 131 — Vl 31 c~ 1 +Vq 31 . All momenta 
are given in GeV. The first entry refers to the standard BT method Eq.(j2U, the 
second to the new method Ea. (j.S2fl and the third to Eq.JSHJ). When no number 
appears in curly brackets, this means that the error does not affect the written 
digits. In the last two cases we have B ~ and in the last case also — m| ~ 0. 
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6. The two-loop vertex V 2S1 

Another imprtant example of two-loop vertex topologies is V 231 (see 
appendix). In this case Feynman parametrisation gives: 



V 



231 



f "-jT(2 + e)JdS 2 (x 1 ,x 2 )[x 2 (l-x 2 )} 
dS 3 (y 1 ,y 2 ,y 3 )y e 3 /2 U~ 2 ~ e , 



-l-e/2 



(37) 



where 



U 



-\P2 2/1 -P(V2- Xy z ) +pi} 2 + (P 2 -pl+mj- m\)y x 
-{P 2 + mg - ml) y 2 + (M 2 - mj) y 3 + pj + m§, 



(38) 



X = 



.(• ] 1 r2 _ -P 2 + ari (P 2 + mf - m|) + x 2 (m 2 - mf) + m| 



x 2 



x 2 (1 - x 2 ) 



(39) 

Now the polynomial U is a quadratic in y±, y 2 and j/3 of the type U = 
y l H y + 2 X* y + L and in principle we could apply the new BT-like relation 
for (i = — 2 — e. Nevertheless, this would not be the most clever approach, 
because in this case the H matrix is singular and therefore the BT-factrors 
B and X are not well-defined. Anyway, thanks to this singularity, we can 
perform the following change of variable 



2/2^2/2 + ^2/3 



(40) 



and obtain a new polynomial U' which is now linear in y 3 . Since this diagram 
is not UV divergent, we can set e = and the integration in y% becomes 
trivial, giving: 



V 



231 




X2 (1 - Xi) 



1 



+ 



L 



dy2 



y 2 )A{x)+x 2 (l-x 1 )B(y) B(y) 
x 2 (xi - x 2 ) 1 



y 2 A{x)+x 2 {x 1 -x 2 )B{y) B{y) 



(41) 
(42) 
(43) 



where A{x) and B{y) are quadratics rispectively in x±, x 2 and y±, y 2 . There- 
fore each term is now a quadratic in y±, y 2 to power —1 with x-dependent 
coefficients. In other words, the y- integrations are one-loop 3-point func- 
tions Co with x-dependent masses and momenta (of course some change of 



10 



variable has to be done to reduce to the usual simplex in y\, y 2 ). Note 
that the zeros of A(x) do not spoil the smoothness of the integrand because 
the C functions cancel one another in the limit A{x) — > 0. At the end the 
diagram is written in the following form 



V 



231 



dS 2 (x 1 ,x 2 } 



1 



A(x) 



r (l) r {2) r (3) Mi) 



(44) 



and computed numerically using for the Co functions the expression of 
Ea. (|24|) obtained with the new method. The numerical results are shown in 
Fig. El 
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Fig. 2. Numerical results for the V 231 family. All momenta are given in GeV. When 
no number appears in curly brackets, this means that the error does not affect the 
written digits. 



7. Work in progress: fermionic correction to sin 2 fllSP * 

The effective leptonic weak mixing angle is at present the observable 
that can give the most stringent indirect evaluation on the mass of the 
Higgs boson (by means of the radiative corrections). In the future, if the 
Higgs boson would be discovered, sin 2 6^ t would represent a strong test of 
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the Standard Model. It is defined throught the vectorial and axial effective 
couplings g v and g a of the Z boson with leptons: 



sin 2 ^ 



4 V \9a 



(45) 



The effective couplings g v and g a are defined throught the amplitude of the 
decay of an on-shell Z boson into two leptons: 



M l zn = u, v, et = u, 7,, (q„ + Qa 7.0 viet, s = Ml 



(46) 



Therefore g v and g a can be obtained from the matrix by using suitable 
projectors: 



1 



1 



TV( 7/ , 7 5M' 



(47) 



where D is the dimension of the space-time. Among the electroweak two- 
loop diagrams contributing to A4^, we started considering those contain- 
ing a closed fermion loop (this computation has been recently done by 
M. Awramik, M. Czakon, A. Freitas and G. Weiglein in ^2])- They are 
represented in Fig. [U 



-y,Z,W 





j,Z,W 




Fig. 3. Diagrams contributing to the fermionic corrections to sin 2 O^S*. Because of 
CP conservation, diagrams involving the Higgs boson cancel and are not included. 

We started appling our methods on configuration c which is the most 
difficult one. The strategy used is the following: 

• Write the amplitude, considering the different contributions to con- 
figuration c. This generates a sum of tensor integral of the family 
V 231 . 

• Perform a simple reduction of tensor integral of the following type: 

2q-p 



{q 2 + m 2 ) [{q + p) 2 + M 2 ] ~ 

1 p 2 — m 2 + M 2 



q 2 +m 2 {q+p) 2 +M 2 (q 2 + m 2 ) [(q + p) 2 + M 2 ] 



(48) 
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This kind of reduction does not introduce any new denominator and 
therefore any spurious singularity. At this point we have a sum of 
scalar and vector integrals belonging to 4 vertex families (V 121 , V 131 , 
V 221 and V 231 ), together with some self-energies and some one- loop 
diagrams. 

• Combine the sum of all these diagrams in just one integral. 

Just to show the efficiency of the numerical computation, we give some 
preliminary results. We consider the expansion in loops for g v /g a 



9v 
9a 



gj 
g° a 



i + 



9v,, 



gl 
g° v 



12 
9v,a 9v,a "T" Qv.a) 



gl 

9°a 




gl 
g° a 



(49) 
(50) 



where the last two terms represent the pure two- loop corrections to sin 2 0Jf 
The contribution to these corrections coming from diagram c (Fig. with 
two Z or two W is: 



(gl 


gl 


u° 


g° a 


gl 


gl\ 


g° v 


g° a 



zz 



WW 



-0.279937 x 10~ 5 ± 0.15 x 10~ 9 , 



0.577269 x 10" 1 ± 0.14 x 10" 



(51) 



In this result are summed up the contributions coming from all possible 
fermion loops. In particular it includes the diagrams containing the top 
quark which in the usual analytical approach require an expansion in the 
ratio Mz/mt ~ 1/4. 



8. Conclusions 

The techniques described in this paper show that the numerical approach 
to Feynman integrals allows the computation of diagrams that cannot be 
treated within the usual analytical approach. Under this point of view, 
would be interesting to apply these methods, and in particular the new one, 
to more complicated two-loop diagrams (the two-loop 4-point functions for 
example). In addition to that, the first results, obtained from the application 
of these techniques to the fermionic correction to sin 2 , seem to show 
that the numerical approach is not only suited for the computation of single 
integrals, but can also be applied to the complete evaluation of physical 
observables (which requires to sum up several diagrams). Of course the 
computaion must be completed to give a serious proof of that. 
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Appendix A 
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Fig. 4. Two-loop vertex topologies. The momenta are all incoming. 



REFERENCES 
[1] [LEP Collaboration] , |arXiv:hep-ex/0312023 

[2] A.I. Davydychev and V.A. Smirnov, Nucl. .lustrum. Meth. A 502 (2003) 621. 
|hep-ph/0210171| . 

[3] R. Bonciani, P. M astrolia and E. Remiddi, Nucl. Phys. B 661 (2003) 289. 
|hep-ph/0301170| . 

[4] F. V. Tkachov, Nucl. Instrum. Meth. A 389 (1997) 309 hep-ph/9609429 ; 
L. N. Bertstein, Functional Analysis and its Applications, 6(1972)66. 

[5] G. Passarino, Nucl. Phys. B 619 (2001) 257 [arXiv:hep-ph/0108252| . 

[6] G. Passarino and S. Uccirati, Nucl. Phys. B 629 (2002) 97 
|arXiv:hep-ph/0112004] . 



14 



[7] A. Ferroglia, M. Passera, G. Passarino and S. Uccirati, Nucl. Phys. B 650 
(2003) 162 |arXiv:hep-p h/0209219 . 

[8] A. Ferroglia, G. Passarino, M. Passera and S. Uccirati, Prepared for 31st 
International Conference on High Energy Physics (ICHEP 2002), Amsterdam, 
The Netherlans; 

A. Ferroglia, G. Passarino, S. Uccirati and M. Passera, Nucl. Instrum. Meth. 
A 502 (2003) 391. 

[9] A. Ferroglia, M. Passera, G. Passarino and S. Uccirati, Nucl. Phys. B 680 



(2004) 199 arXiv:hep-ph/0311186 



[10] S. Actis, A. Ferroglia, G. Passarino, M. Passera and S. Uccirati, 
|arXiv:hep~ph/0402T32| 

[11] R. Erdelyi et al., Higher Transcendental Functions vol. 2, Bateman Manuscript 
Project (McGraw-Hill 1953); 

L. J. Slater, Generalized Hypergeometric Functions Cambridge Univ. Press, 
1966. 



[12] M. Awramik, M. Czakon, A. Freitas and G. Weiglein, arXiv:hep-ph/0407317 



